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INTRODUCTION 


An  earlier  report1  attempts  to  reconcile  the  multipolar  lattice 
sums,  Anm,  of  various  structures  to  the  phenomenological  crystal  field 
parameters,  Bnm,  for  a  number  of  rare-earth  ions.  In  the  theory  that 
was  used,  three  parameters  were  introduced.  Par  calcium  tungstate,  the 
three  parameters  introduced  were  t,  the  outward  expansion  of  the  radial 
wave  functions  from  the  Hartree-Fock  value;  q,  the  effective  charge  on 
the  oxygen  ion;  and  g,  an  effective  position  of  the  oxygen  with  respect 
to  the  tungsten  site  chosen  such  as  to  reproduce  the  lowest  of  the 
(W04)-2  complex  when  viewed  from  a  distance  large  compared  with  the  W-0 
distance.  This  model,  when  used  with  appropriate  shielding  factor^,2-6 
gave  excellent  results  for  calcium  tungstate.  Later  improvements6  in 
the  three-parameter  theory  gave  better  results  for  calcium  tungstate. 
Fbr  other  host  materials  such  as  YVOj,  the  same  procedure  does  not  give 
satisfactory  results. 

Therefore,  it  was  decided  to  investigate  the  possibility  of  a 
correction  due  to  the  penetration  of  the  rare-earth  electrons  into  the 
charge  distribution  on  the  ligands.  The  exchange  of  rare-earth  elec¬ 
trons  with  the  electrons  on  the  ligands  was  ignored  in  the  investi- 


1 Richard  P.  Leavitt,  Clyde  A.  Morrison,  and  Donald  E.  Wortman,  Rare 
Earth  Ion-Host  Crystal  Interactions,  3,  Three-Parameter  Theory  of 
Crystal  Fields,  Harry  Diamond  Laboratories  HDL-TR-1673  (June  1975). 

2R.  M.  Sternheimer ,  Phys.  Rev.,  84_  (1951),  244. 

3R.  E.  Watson  and  A.  J.  Freeman,  Phys.  Rev.,  135  (1964),  A1209. 
kD.  Sengupta  and  J.  o.  Artman,  Phys.  Rev.,  Bl_  (1970),  2986. 

5P.  Erdos  and  J.  H.  Kang,  Phys.  Rev.,  B6_  (1972),  3393. 

6 Clyde  A.  Morrison,  Nick  Karayianis,  and  Donald  E.  Wortman,  Rare-Earth 
Ion-Host  Lattice  Interactions,  4.  Predicting  Spectra  and  Intensities  of 
Lanthanides  in  Crystals,  Harry  Diamond  Laboratories  HDL-TR-1816  (June 
1977). 
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gation.  The  effects  of  charge  penetration  were  previously  calculated, 
but  the  method  followed  did  not  leave  any  parameters  that  could  be  used 
to  fit  the  experimental  data?  thus,  this  method  is  not  suitable  for  our 
purposes . 

2 .  THEORY 

The  charge  distribution,  p(r^),  surrounding  the  ligand  ion  is  chosen 
to  be  spherically  symmetric  and  of  the  form 

-Yr, 

P^r^J  =  Nr^e  ,  ( 1 ) 

where  N  =  y5Q/96n  and  Q  is  the  effective  charge  on  the  ligand.  Such  a 
charge  density  is  expected  for  hydrogenic  wave  functions  representing 
the  six  2p  electrons  or  for  both  the  2s  and  2p  electrons  using  five 
Slater-type  orbitals.8 

We  let  R  be  the  position  of  the  nucleus  of  a  ligand  ion  measured 
from  the  rare-earth  nucleus  and  r  be  the  position  of  the  rare-earth 
electron.  Also,  we  let  r^  be  the  position  of  the  volume  element  of 
charge  p(r.j)  dt1  from  the  oxygen  nucleus.  Then  the  electric  potential 
at  the  rare-earth  electron  at  r  due  to  this  element  of  charge  is 


where 


d$ 


111 


**  -► 
Ri  +  rt 


(2) 


?A.  K.  Raychaudhuri  and  D.  K •  Ray,  Proc.  Phys .  Soc.,  90_  (1967),  839, 
eJ.  C.  Slater,  Phys,  Rev,,  36_  (1930),  57. 
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1 


+  ♦  ♦ 

R1  =  R  -  r  , 

<3t^  =  dr^  dft.j  , 

and 

dfi^  =  sin  6^  d0^  d<{) ^  . 


Before  integrating  equation  (2),  it  is  convenient  to  expand  |  +  r^  in 

spherical  harmonics  so  that 


4it 


R1  +  r 


-¥ 

r1 


Jim 


2JI  +  1 


(-1) 


Ml 


Jim 


,3> 


for  r1  <  R^  j  for  r1  >  R^,.  the  factor  in  brackets  becomes  R^/r^+1.  Since 
we  assumed  that  the  charge  density  p(rj)  is  spherically  symmetric,  we 
can  integrate  over  dfi.,  so  that 


for  R1  >  r1 
for  R1  <  r1 


(4) 


where  we  have  used  the  relation 
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The  electric  potential  at  the  rare-earth  electron  is  then 


The  first  integral  in  equation  (6)  gives  the  total  charge  and,  since 
R.j  <  r^,  the  integrand  of  the  second  integral  is  always  positive.  Thus, 
the  potential  is  always  less  than  that  of  a  point  charge  located  at  . 
If  we  let 


I  = 


(7) 


and  use  the  charge  distribution  given  in  equation  (1),  this  becomes 


N4 1  l  3  Y?R1  y3R^ 

I=37l1  +  4YRi  +  — +  -irj 


-YR, 


(8) 


where  we  have  used  the  easily  derivable  expression 


r  n  ~Yx 

I  x  e  dx 

J  y 


nle 


-yy 


n+1 


Y 


8 


-YR« 

Further,  by  differentiating  e  with  respect  to  y,  the  result  on  the 
right  side  of  equation  (8)  can  be  written 


+ 


Y2R2 

4 


+ 


Y3r3\  -yRi 

— r 


3  d  y2 

=  1  —  —  Y  -  +  J — 

4  '  dy 


d2  y 3  <j3  e 


-YR. 


dy 


y  dJ 

2  "  24  dy3 


■YR, 


=  0(Y) 


(9) 


and  equation  (8)  becomes 


I 


-YR, 

e 


The  entire  expression  of  equation  (6)  can  be  written 


4> 


(10) 


where  R,  =  |  R  -  r  | . 

It  is  convenient  for  later  work  to  expand  equation  (10)  in  terms  of 

A 

Cnm*r^'  defined  by 


SunW  '  (^5r)l/2 


9 


9 


so  that  the  first  term  in  equation  (10)  becomes 


n,m 


(11) 


with  r  <  R.  Similarly,5 


(2n+1)i  ( yr  )k  (yR)C*  (r)c  (r) 
7T  n  n  nmx  '  nm'  ' 

n,m 


(  12) 


where  R  <  r  and  in(x)  and  kn(x)  are  modified  spherical  Bessel  functions 
of  the  first  and  third  kinds,  respectively. 


If  we  define 


A 


nm 


and 


A  (r,R)  =  - -  2(2naM)  0(y)yi  { yr  )k  (yR) 

n  n  Ti  n  n 

r 

we  can  write  the  potential  energy  of  the  electron  at 
(10)  as 


,  (13) 

r  using  equation 


U(r)  ■  -e<)>(r) 


(14) 


EA*  [l  -  A  (r,R)l  rnC  (r) 
nm  I  n  I  nmv  ' 

n  .in  u  J 


n#m 


Abramcwitz  and  I.  A.  Stegun,  Handbook  of 
National  Bureau  of  Standards  Applied  Mathematics 


Mathematical  Functions , 
Series  (1966),  445. 
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which  shows  that,  if  An(r,R)  is  positive  and  less  than  1,  the  charge 
distribution  reduces  the  Anm  that  would  be  calculated  from  a  point 
charge  model.  Since,  as  R  becomes  very  large,  kn(yR)  decreases  as  e-^R, 
An(t,R)  vanishes,  and  the  rare-earth  electron  "sees"  effectively  a  point 
charge.  Thus,  for  ions  remote  from  the  rare-earth  site,  a  point  charge 
model  is  sufficient. 


3.  DISCUSSION  OF  AND  APPROXIMATIONS  TO  An(r,R) 

In  defining  An(r,R)  by  equation  (13),  the  expression 

°(Y)  [\i.n<  Yr)kn(  yR)  J 


occurs.  With  a  great  deal  of  tedious  algebra  using  O(y)  defined  by 
equation  (9)  and  using  the  differential  equations  satisfied  by  in(x)  and 
kn(x),  this  result  is  given  by 


0(y)yi  <yr)k  (yR)  =  (2n+1^ -  f F  ( x ,  y )  i  (x)k  (y) 

n  n  „  n  in  nn 

2irx 


-  Gn(x'y)in-1(x>kn(y)  +  Vx'y)in(x)kn-1(y) 


(15) 


-3xyin.1(X,kn-1(y)]  ' 


with  x  =  yr  and  y  =  yR.  In  this  derivation,  the  following  were  used: 
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r 


i 


and 


z2w"  +  2zw'  -  fz2  +  n(n+1)]w  =  0  , 

n  n  L  ‘  n 


with  wn  either  in  or  kn  and  z  given  by  x  or  y,  respectively.  The 
various  functions  in  equation  (15)  are  defined  by 


F  (x,y )  =  (2n+3)(nt2)(2n+5)  +  4n+9  (x2+y2; 


.  .  .  ( 2n2+1 ln-t-18  .  y2  .  x2 \ 

;n(«.y)  - - - **r  +  T)  ■ 


H  (x,y)  =  y 
n 


^2n2+1 


^1-18  +  >c£  + 


2  6 


(  16) 


The  result  given  in  equation  (15),  which  is  rigorously  correct,  is  not 
very  desirable.  In  the  theory  of  crystal  fields  using  the  point  charge 
model,  the  Anm  are  used  in  conjunction  with  shielding  factors  and 
appropriate  radial  wave  functions  to  attempt  to  fit  the  phenomenological 
crystal  field  parameters,  Bnm,  which  are  obtained  by  fitting  the  crystal 
field  spectra  of  a  given  ion  by  using  the  hamiltonian 


H  ~  5""'  B*  C  (r)  •  (  17) 

i—t  nm  nm  '  ' 
n,m 

With  the  inclusion  of  the  factor  An(x,y)  given  in  the  form  of  equation 
(15),  this  combination  would  be  possible,  but  would  be  extremely 
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cumbersome.  Fbr  this  reason,  we  have  chosen  to  replace  An(x,y)  by  an 
approximation  for  small  x  (small  yr). 


Before  proceeding,  it  is  convenient  to  make  a  few  simplifying 
changes  in  variables.  We  let 

.  2nxnn! 

1n  X  ( 2n+ 1 ) !  Rn  ' 


.  ,  ,  Tte  y  ( 2n)  ! 

k  (y)  =  - 7  -  q 

n  „  n+1  n  ,  n 

2y  2  n ! 


where 


(  18) 


Rn  1  +  2 ( 2n+3 )  +  8 ( 2n+3 ) < 2n+5 ) 


+  .  . 


<30  -  1 


q,  =  1  +  y 


%  =  V 


+  - 


y2V2 


1  (2n-1)(2n-3) 


With  these  changes,  An(x,y)  (eq  15)  becomes 
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(19) 


n(x,y)  =  ^[Fn(x,y)qnRn  -  Gn(x,y)  ^  q^., 


+  Hn(X'y)  2^  VlRn 


->  2  2n+1  1 

y  2n-1  qn-1  n-1 J 


In  this  form,  it  is  apparent  from  equations  (16)  and  (19)  that  An(x,y) 
is  a  function  of  x2  and,  further,  the  limit  x  =  0,  y  =  0  gives 
An(0,0)  =  1.  For  small  x,  the  following  expansion  for  An(x,y)  can  be 
obtained : 


A  ( x ,  y ) 
n 


A  ( 0 , y )  +  x2A* (0,y)  +  x4A"(0,y) 
n  n  n 


(20) 


where 


A  (0,y) 
n 


j[i  +  T  y2]  qn  +  (2n2 '  7n  +  9  +  rk-l)  ' 


2  -y  r  2  ”1 

*;<°'y)  =  24(2n«)  “nqn  +  [n  +  ' 


An(0,y>  96  (  2n+3 )  T 2n+5 )  i-(n+3)c* 


{-(n+3)qn  +  2n-T  _( 


(n+1 ) ( 2n+3 )  + 


}%-]  ■ 


The  first  term  in  An(x,y)  from  equation  (20)  is,  in  principle,  only  host 
dependent  as  y  =  yR  and  y  is  dependent  on  the  ligand  ion,  whereas  R  is 
the  distance,  which  in  the  absence  of  local  distortion  is  characteristic 
of  the  host  material.  The  other  two  terms  in  equation  (20)  are 
dependent  on  the  particular  rare-earth  ion  through  the  factors  x2  and 
x*4.  From  equation  (14),  we  see  that  the  total  radial  dependence  of  the 
term  containing  An(x,y)  is  of  the  form  rnAn(x,y)  and,  with  x  =  yr,  for 
the  equivalent  electrons,  we  can  write 
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<rn&  (x,y)>  =  <r°>A  (0,y)  +  Y2<rn+2>A' ( 0,y ) 
n  n  n 


(21) 


+  Y<4<rn+4>A,,(0,y) 
n 


k  k 

where  <r  >  is  the  expectation  value  of  r  for  the  rare-earth 

electrons.  Frequently,  Hartree-Fock  wave  functions  modified  in  some 

manner  are  used  to  calculate  these  expectation  values.  However,  the 

unmodified  Hartree-Bock  functions  may  be  sufficiently  accurate  to  allow 

us  to  write  equation  (21)  as 


* 


<rnA  (x,y)>  =  <rn>  [a  (0,y)  +  Y2X_(n)A* (0,y)  +  Y4*,  (n) A"( 0,y)l  ,  (22) 

n  Ln  z  n  4  n  j 


where 


X2(n)  =  <rn+2>/<rn>  , 


A  (n)  =  <rn+4>/<rn> 

4 


with  the  A2(n)  aid  A^(n)  calculated  once  and  for  all  for  each  rare-earth 
ion  by  using  Hartree-Fock  wave  functions.  The  modifications  of  these 
functions  could  then  be  employed  in  the  <rn>.  The  functions  A2(n)  and 
A4(n)  were  calculated  for  a  number  of  rare-earth  ions,  and  the  results 
for  Nd3+  and  Er3+  are  given  in  appendix  A. 
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4.  ESTIMATION  OF  y  and  R 

In  theory,  the  functions  An(0,y),  A^(0,y),  and  A^(0,y)  are  dependent 
only  on  the  host  lattice  through  y,  which  describes  the  size  of  the 
ligand  charge  distribution,  and  the  distance,  R,  from  the  rare-earth 
site  to  the  nearest  ligand.  It  would  be  advantageous  to  have  these 
quantities  tabulated  for  a  range  of  values  of  y  and  R  centered  at  some 
typical  values.  A  rough  estimate  of  y  for  oxygen  can  be  obtained  by 
using  the  charge  distribution  given  by  equation  ( 1 )  to  calculate  the 
average  value  of  r  by 

Jrp  dx  5 
r  =  — c -  =  —  , 

Q  y 

and  if  we  let  r  be  the  ionic  radius  of  oxygen,  which  is  approximately 
140  nm,  we  obtain 


y  ~  3 • 57/A  . 

A  second  estimate  of  y  for  oxygen  can  be  obtained  by  using  rules  given 
by  Slater8  for  'deriving  the  exponential  factors  in  the  wave  functions 
for  the  2p8  (also  the  2s2)  of  oxygen  to  obtain  a  value 

y  ~  3.86/A  , 

where  this  value  of  y  is  twice  that  given  by  Slater's  rule  since  the 
charge  distribution  is  proportional  to  the  square  of  the  wave 
function.  To  obtain  an  estimate  of  typical  values  of  R,  the  x-ray  data 


8J.  C.  Slater ,  Phys.  Rev.,  36  (1930),  57. 
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on  a  number  of  crystals  were  examined.  The  smallest  distance  found  was 

R^2.1  A  (the  smallest  Yb-0  distance  in  YbPjOg),  and  the  largest 

distance  found  was  near  R  =  2.6  A  (the  smallest  Ho-0  distance  in 

HoV04).  A  further  check  to  find  a  typical  range  for  y  =  yR  was  an 

approximate  fitting  of  the  theory  to  the  phenomenological  Bnm,  which 

gave  y  ~  10.  It  was  then  decided  that  the  functions  ^(O.y), 

A'(0,y),  and  A"(0,y)  should  be  calculated  for  the  range  of  values 

n  n 

8  <  y  <  10  (fig.  1  to  3).  The  results  for  y  =  9.6,  10,  and  10.4  are 

given  in  tabular  form  in  appendix  B.  Several  features  of  these 

functions  are  of  interest.  The  function  An(0,y)  (fig.  1)  always 

increases  for  larger  n  with  y  fixed,  whereas  the  function  A* (0,y)  (fig. 

n 

2)  is  larger  for  intermediate  n  (n  ~  4)  for  small  y,  changing  when  y 

gets  larger.  The  function  A"(0,y)  (fig.  3)  is  negative  for  some  values 

n 

of  n  for  y  <  10.4  and  is  positive  for  larger  y.  Each  correction  is  of 
the  order  of  10~2  of  the  preceding  term;  for  example. 


A^(0,y)  ~  An(0,y)  x  10~2 


4, 
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5,  ,  CONCLUSION 


u 


Y 

l/t 


The  theory  of  charge  perpetration  by  a  rare-Mrth  ion  into  the  charge 
distribution  on  a  ligand  )tas  been  developed,  h he  results  are  cast  into 
a  form  reminiscent  of  the  Sternheimer  shielding  factors,  which 

,  n.--  /  r)  /  n 

relate  the  observed  Ay7<r  >  to  the  lattice/sum  A__<r  >  in  the  form 

nm  *  /■'  ,  i 

i  A  .  >  ^ 

*  ^  ^ ; 


A^i,<rn>  =  A^<rn>/1  - 
nm  r>ro  \  n 


(23) 


The  corresponding  form  developed  here  is 


/ 


) 


v  *  > 


-  »»«"> [’ - #»•*>/«”>]  •  ' 


(24) 


with  explicit  expressions  given  for  the  A^(x,y).  The  zi^lx,y)  is  then 
expanded  for  small  x,  and  a  form  convenient  for  calculation  is  given. 


The  results  obtained  here  can  be  used  with  the  three-parameter 

'-~7 

1  a —  — %_j —  .i  — - -  i  Ci  j  _  -  — - ~  —  n  if 


theory  of  crystal  fields1  to  obtain  crystal  field  parameters  Bnnl(  Tjqny)  # 
which  give  a  better  representation  of  the  phenomenological  t^ian  i\e 
obtained  previously. 
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1 Richard  P.  Leavitt,  Clyde  A.  Morrison,  and  Donald  E.  Wortman,  Pare 
Earth  Ion-Host  Crystal  Interactions,  3.  Three-Parameter  Theory  of 
Crystal  Fields,  Harry  Diamond  Laboratories  HDL-TR-167 3  (June  1975). 
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APPENDIX  A 


In  equation  (22)  in  the  main  body  of  this  report,  the  quantities 
X2(n)  =  <rn+2>/<rn>  and  X4(n)  =  <rn+4>/<rn>  are  introduced.  These 

quantities  are  given  in  tables  A-1  and  A-2  for  Nd3+  and  Er3+, 
respectively.  The  Hartree-Fbck  wave  functions  of  Freeman  and  Watson* 
are  used. 


TABLE  A-1.  VALUES  OF  X2(n)  =  <rn+2>/<rn>  AND  X4(n)  =  <rn+4>/<rn> 
FOR  Nd3+ 


n 

X2(n) 

(A) 

X4<n) 

(A) 

0 

0.28030 

0.18831 

1 

0.43768 

0.43987 

2 

0.67181 

0.97107 

3 

1.00500 

1.98782 

4 

1.44545 

3.71583 

5 

1.97793 

6.31501 

6 

2.57070 

9.83261 

7 

3.19274 

14.23983 

8 

3.82487 

19.50188 

TABLE  A-2.  VALUES  OF  X2(n)  =  <rn+2>/<rn>  AND  X4(n)  =  <rn+4>/<rn> 
FOR  Er3+ 


n 

X2(n) 

(A) 

A4Cn) 

(A) 

0 

0.18652 

0.08828 

1 

0.30141 

0.21297 

2 

0.47331 

0.46836 

3 

0.70659 

0.91977 

4 

0.98955 

1.60912 

5 

1.30171 

2.54414 

6 

1.62611 

3.71580 

7 

1.95446 

5.11960 

8 

2.28508 

6.76357 

*A.  J.  Freeman  and  R.  E.  Watson,  Phys. 


Rev,, 


127  (1962),  2058. 
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APPENDIX  B 


In  equations  (22)  and  (23)  in  the  main  body  of  this  report,  the 
functions  An(0,y),  A^(0,y),  and  AB(0,y)  are  introduced  and  depend  on  the 
spread  of  the  wave  function,  y,  of  the  ligand  and  the  distance,  R,  to 
the  ligand  through  y  =  yR.  In  the  main  body,  it  is  shown  that  a  value 
of  y  ~  10  is  representative  for  oxide  compounds.  These  three  quantities 
are  given  for  y  =  9.6,  10.0,  and  10.4  in  tables  B-1,  B-2,  and  B-3, 
respectively. 


TABLE  B-1.  FUNCTIONS  A„(0,y),  A'(0,y),  AND  AB(0,y)  PDR  y  =  9.6 

n  n  n 


n 

An(0,y) 

A^  (  0  *  y ) 

An  ( o ,  y ) 

i 

3.7795(-2) 

1.8975(-3) 

2.5413(-5) 

2 

1. 1450(-1) 

3. 7665 (-3) 

3. 1391 ( -5) 

3 

2.31 08 ( — 1 ) 

5.3S84(-3) 

2.6028(-5 ) 

4 

3.6993(-1 ) 

6.2777(-3 ) 

1.3800 (-5) 

5 

5. 1102 (-1) 

6.4717(-3) 

-5.3134(-8) 

6 

6.3959(-1 ) 

6.0987 (-3 ) 

-1. 1995 (-5) 

Note:  Numbers  in  parentheses  are  powers 
of  10. 


TABLE  B-2.  FUNCTIONS  An(0,y),  A^(0,y),  AND  A£(0,y)  FOR  y  =  10 


n 

An(0,y) 

A^(0,y) 

A£(0,y) 

i 

2.9253(-2) 

1 . 5 133 (-3 ) 

2. 1619(-5) 

2 

9.2308(-2 ) 

3. 1528 (-3 ) 

2.9026(-5) 

3 

1.9320(-1) 

2.7006(-5) 

2.7006(-5) 

4 

3. 1931 (-1 ) 

1.8079(-5) 

1.8079(-5) 

5 

4.5343(-1 ) 

6.3338(-6) 

6.3338(-6) 

6 

5.8099(-1 ) 

-4.9131 (-6) 

-4.9131 (-6) 

Note:  Numbers  in  parentheses  are  powers 
of  10. 
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TABLE  B-3.  FUNCTIONS  A„(0,y),  A'(0,y),  AND  A"(0,y)  FOR  y  =  10.4 

11  n  n 


n 

V0,y) 

A^(0,y) 

A" ( 0 ,y ) 
n 

i 

2.2532(-2) 

1.1981 (-3) 

1 .8 1 04 ( -5  ) 

2 

7.3957 ( -2 ) 

2.6136(-3) 

2.6144( -5 ) 

3 

t .6035 (-1 ) 

4 . 0634 ( -3  ) 

2 .6664 (-5 ) 

4 

2. 7343 (-1  ) 

2.0692(-5) 

5.  1755(-3 ) 

5 

3.9900(-1  ) 

5.7737 (-3 ) 

1 . 1 199 ( -5 ) 

6 

5.2339{-1 ) 

1.0883 (-6) 

5.8660 ( -3 ) 

Note:  Numbers  in  parentheses  are  powers 
of  10. 
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